Role of SARS-CoV-2 mutations in the evolution of the COVID-19 pandemic.

OBJECTIVES: The SARS-CoV-2 pandemic and large-scale genomic surveillance provided an exceptional opportunity to analyze mutations that appeared over three years in viral genomes. Here we studied mutations and their epidemic consequences for SARS-CoV-2 genomes from our center. METHODS: We analyzed 61,397 SARS-CoV-2 genomes we sequenced from respiratory samples for genomic surveillance. Mutations frequencies were calculated using Nextclade, Microsoft Excel, and an in-house Python script. RESULTS: A total of 22,225 nucleotide mutations were identified, 220 (1.0%) being each at the root of ≥836 genomes, classifying mutations as 'hyperfertile'. Two seeded the European pandemic: P323L in RNA polymerase, associated with an increased mutation rate, and D614G in spike that improved fitness. Most 'hyperfertile' mutations occurred in areas not predicted with increased virulence. Their mean number was 8±6 (0-22) per 1000 nucleotides per gene. They were 3.7-times more frequent in accessory than informational genes (13.8 versus 3.7/1000 nucleotides). Particularly, they were 4.1-times more frequent in ORF8 than in the RNA polymerase gene. Interestingly, stop codons were present in 97 positions, almost only in accessory genes, including ORF8 (21/100 codons). CONCLUSIONS: most 'hyperfertile' mutations did not predict emergence of a new epidemic, and some were stop codons indicating the existence of so-named 'non-virulence' genes.

Vγ9Vδ2 T-cells are potent inhibitors of SARS-CoV-2 replication and represent effector phenotypes in COVID-19 patients.

Vγ9Vδ2 T-cells play a key role in the innate immune response to viral infections through butyrophilin (BTN)-3A. Here, we reported that blood Vγ9Vδ2 T-cells decreased in clinically mild COVID-19 compared to healthy volunteers (HV), and was maintained up to 28-days and in the recovery period. Terminally differentiated Vγ9Vδ2 T-cells tend to be enriched on the day of diagnosis, 28-days after and during the recovery period. These cells showed cytotoxic and inflammatory activities following anti-BTN3A activation. BTN3A upregulation and Vγ9Vδ2 T-cell infiltration were observed in a lung biopsy from a fatal SARS-CoV-2 infection. In vitro, SARS-CoV-2 infection increased BTN3A expression in macrophages and lung cells that enhanced the anti-SARS-CoV-2 Vγ9Vδ2 T-cells cytotoxicity and IFNγ and TNFα. Increasing concentrations of anti-BTN3A lead to viral replication inhibition. Altogether, we report that Vγ9Vδ2 T-cells are important in the immune response against SARS-CoV-2 infection and that activation by an anti-BTN3A antibody may enhance their response.

Incidental diagnosis of mpox virus infection in patients undergoing sexually transmitted infection screening-findings from a study in France.

OBJECTIVES: This study aimed to investigate the prevalence of mpox virus (MPXV) infections in the general population consulting for routine sexually transmitted infections (STIs) in our Marseille public hospital. In fact, the recent worldwide MPXV outbreak mainly impacted men who have sex with men and the prevalence of MPXV infections in the general population remains poorly defined. METHODS: All samples addressed routinely to our microbiological laboratory for STIs between July 1 and October 15, 2022 were screened with MPXV-specific quantitative polymerase chain reaction. RESULTS: A total of 2688 samples from 1896 patients suspected of having STIs were tested and eight (0.4%) patients were incidentally diagnosed with MPXV infection, including six men and two women. MPXV was detected in rectal swabs (n = 2), urine (n = 2), vaginal swabs (n = 2), a urethral swab (n = 1), and a skin swab (n = 1). CONCLUSIONS: This study suggests that some MPXV infections are likely to be underdiagnosed because of their non-specific clinical presentation and/or insufficient clinical knowledge of the disease. Our data showed that systematic screening was particularly useful for detecting MPXV in patients without classic lesions or cases of asymptomatic carriage in patients reporting recent high-risk exposure and in patients presenting no obvious risk factor.

Resistance-associated mutations to the anti-SARS-CoV-2 agent nirmatrelvir: Selection not induction.

Mutations associated with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) resistance to antiprotease nirmatrelvir were reported. We aimed to detect them in SARS-CoV-2 genomes and quasispecies retrieved in our institute before drug availability in January 2022 and to analyze the impact of mutations on protease (3CLpro) structure. We sought for 38 3CLpro nirmatrelvir resistance mutations in a set of 62 673 SARS-CoV-2 genomes obtained in our institute from respiratory samples collected between 2020 and 2023 and for these mutations in SARS-CoV-2 quasispecies for 90 samples collected in 2020, using Python. SARS-CoV-2 protease with major mutation E166V was generated with Swiss Pdb Viewer and Molegro Molecular Viewer. We detected 22 (58%) of the resistance-associated mutations in 417 (0.67%) of the genomes analyzed; 325 (78%) of these genomes had been obtained from samples collected in 2020-2021. APOBEC signatures were found for 12/22 mutations. We also detected among viral quasispecies from 90 samples some minority reads harboring any of 15 nirmatrelvir resistance mutations, including E166V. Also, we predicted that E166V has a very limited effect on 3CLpro structure but may prevent drug attachment. Thus, we evidenced that mutations associated with nirmatrelvir resistance pre-existed in SARS-CoV-2 before drug availability. These findings further warrant SARS-CoV-2 genomic surveillance and SARS-CoV-2 quasispecies characterization.

The return of the "Mistigri" (virus adaptative gain by gene loss) through the SARS-CoV-2 XBB.1.5 chimera that predominated in 2023.

Severe acute respiratory syndrome coronavirus 2 XBB.1.5 is the first recombinant lineage to predominate at the country and global scales. Very interestingly, like the Marseille-4B subvariant (or B.1.160) and the pandemic variant B.1.1.7 (or Alpha) previously, it has its ORF8 gene inactivated by a stop codon. We aimed here to study the distribution of stop codons in ORF8 of XBB.1.5 and non-XBB.1.5 genomes. We identified that a stop codon was present at 89 (74%) ORF8 codons in ≥1 of 15 222 404 genomes available in GISAID. The mean proportion of genomes with a stop codon per codon was 0.11% (range, 0%-7.8%). In addition, a stop codon was detected at 15 (12%) codons in at least 1000 genomes. These 15 codons are notably located on seven stem-loop hairpin regions and in the signal peptide region for the case of the XBB.1.5 lineage (codon 8). Thus, it is very likely that stop codons in ORF8 gene contributed on at least three occasions and independently during the pandemic to the evolutionary success of a lineage that became transiently predominant. Such association of gene loss with evolutionary success, which suits the recently described Mistigri rule, is an important biological phenomenon very unknown in virology while largely described in cellular organisms.

How to distinguish mpox from its mimickers: An observational retrospective cohort study.

During the current global outbreak of mpox (formerly monkeypox), atypical features were frequently described outside endemic areas, raising concerns around differential diagnosis. In this study, we included 372 adult patients who had clinical signs consistent with mpox and who were screened using non-variola orthopoxvirus specific quantitative polymerase chain reaction (PCR) between 15 May and 15 November 2022 at the University Hospital Institute Méditerranée Infection, Marseille, France. At least one clinical sample was positive for 143 (38.4%) of these patients and 229 (61.6%) were negative. Clinically, patients who had mpox presented more frequently with systemic signs (69.9% vs. 31.0%, p < 10-6 ) including fever (51.0% vs. 30.1%, p < 10-3 ), myalgia (33.5% vs. 17.9%, p = 0.002), and lymphadenopathy (38.5% vs. 13.1%, p < 10-6 ). Among the patients who were negative for the non-variola orthopoxvirus, an alternative diagnosis was identified in 58 of them (25.3%), including chickenpox (n = 30, 13.1%), syphilis (n = 9, 4%), bacterial skin infection (n = 8, 3.5%), gonococcus (n = 5, 2.2%), HSV infection (n = 5, 2.2%), and histoplasmosis (n = 1, 0.4%). Overall, in the current outbreak, we show that mpox has a poorly specific clinical presentation. This reinforces the importance of microbiological confirmation. In symptomatic patients who are negative for the monkeypox virus by PCR, a broad differential diagnosis should be maintained.

Leucobacter manosquensis sp. nov.-A Novel Bacterial Species Isolated from Healthy Human Skin.

Extending our knowledge on human skin microbiota is a challenge to better decipher its role in health and disease. Using the culturomics method, we isolated strain Marseille-Q4368 from the healthy forehead of a 59-year-old woman. We describe here the main characteristics of this bacterium using a taxonogenomic approach. This new bacterial species is Gram-positive, non-motile, and non-spore-forming. Its 16S rRNA sequence exhibited a similarity of 99.59% with Leucobacter chromiiresistens, the most closely related species in terms of nomenclature. However, a digital DNA-DNA hybridization analysis between these two species revealed a maximum identity similarity of only 27.5%. We found phenotypical and genomic differences between strain Marseille-Q4368 and its closely related species. These findings underscore the classification of this bacterium as a distinct species. Hence, we propose the name Leucobacter manosquensis sp. nov. strain Marseille-Q4368 (=CSUR Q4368 = DSM 112403) for this newly identified bacterial species.

Emergence of a second SARS-CoV-2 variant with a tremendous genetic leap from its ancestors.

The on-going emergence of the Omicron BA.2.86 variant is one of the major events in SARS-CoV-2 genetic evolution that remain enigmatic regarding the overall virus mutation rate, together with the emergence of the initial Omicron variant, BA.1. Indeed, the genomes of the Omicron BA.2.86 lineage, an offspring of the second major Omicron subvariant, BA.2, harbor 39 additional mutations in the spike compared to this ancestor. Here we comment on the phylogeny of BA.2.86, on the positions, and frequencies in other SARS-CoV-2, of mutations in its spike, and on the structural model of this protein that concentrates most of BA.2.86 additional mutations.

Temporal Shift of the Respiratory Syncytial Virus Epidemic Peak for the Years 2020-2023 in Marseille, Southern France.

Respiratory syncytial virus is among the most common causes of respiratory infections. Typically, this viral infection has a seasonality during the cold months but with the SARS-CoV-2 pandemic this has been considerably modified. Here, we studied the epidemiology of this virus in university hospitals of Marseille, South of France, over the period 2020 to 2023. We tested in our laboratory from July 2020 to October 2021 16,516 nasopharyngeal swabs from 16,468 patients for RSV infection using different qPCR assays. We then analyzed data from previous and subsequent winters (from 2018 to 2023) and previous summers (from 2015 to 2021). A total of 676 patients were RSV-positive; their mean age was 3 years and 91 were under 5 years of age. We observed a delay of 4 months of the RSV epidemic's onset compared to other years with an epidemic that peaked in March 2021. We had significantly more RSV-positive cases during summer 2021 compared to previous summers, whereas the incidence of RSV infections was not significantly higher during winter 2022 versus previous winters. Moreover, 494 patients were diagnosed as RSV-positive in the emergency unit and 181 were subsequently hospitalized, and 34 patients were diagnosed RSV-positive while already in the intensive care unit. Over all the study periods, 38 patients diagnosed as RSV-positive died, the majority of whom (23/28) were over 65 years of age. These data show an atypical evolution of the incidence of RSV infections in our city and is another example of the unpredictability of infectious disease epidemiology.

New putative periodontopathogens and periodontal health-associated species: A systematic review and meta-analysis.

To investigate the existence of any association between new putative periodontal pathogens and periodontitis. Two independent reviewers conducted electronic literature searches in the MEDLINE (PubMed), EMBASE, DOSS and Google Scholar databases as well as a manual search to identify eligible clinical studies prior to November 2022. Studies comparing the prevalence of microorganisms other than the already-known periodontal pathogens in subgingival plaque and/or saliva samples between subjects with periodontitis and subject with periodontal health were included. Meta-analyses were performed on data provided by the included studies. Fifty studies including a total of 2739 periodontitis subjects and 1747 subjects with periodontal health were included. The Archaea domain and 25 bacterial species (Anaeroglobus geminatus, Bacteroidales [G-2] bacterium HMT 274, Desulfobulbus sp. HMT 041, Dialister invisus, Dialister pneumosintes, Eubacterium brachy, Enterococcus faecalis, Eubacterium nodatum, Eubacterium saphenum, Filifactor alocis, Fretibacterium sp. HMT 360, Fretibacterium sp. HMT 362, Mogibacterium timidum, Peptoniphilaceae sp. HMT 113, Peptostreptococcus stomatis, Porphyromonas endodontalis, Slackia exigua, Streptococcus gordonii, Selenomonas sputigena, Treponema amylovorum, Treponema lecithinolyticum, Treponema maltophilum, Treponema medium, Treponema parvum and Treponema socranskii) were found to be statistically significantly associated with periodontitis. Network studies should be conducted to investigate the role of these newly identified periodontitis-associated microorganisms through interspecies interaction and host-microbe crosstalk analyses.

Sequencing of monkeypox virus from infected patients reveals viral genomes with APOBEC3-like editing, gene inactivation, and bacterial agents of skin superinfection.

A large outbreak of Monkeypox virus (MPXV) infections has arisen in May 2022 in nonendemic countries. Here, we performed DNA metagenomics using next-generation sequencing with Illumina or Nanopore technologies for clinical samples from MPXV-infected patients diagnosed between June and July 2022. Classification of the MPXV genomes and determination of their mutational patterns were performed using Nextclade. Twenty-five samples from 25 patients were studied. A MPXV genome was obtained for 18 patients, essentially from skin lesions and rectal swabbing. All 18 genomes were classified in clade IIb, lineage B.1, and we identified four B.1 sublineages (B.1.1, B.1.10, B.1.12, B.1.14). We detected a high number of mutations (range, 64-73) relatively to a 2018 Nigerian genome (genome GenBank Accession no. NC_063383.1), which were harbored by a large part of a set of 3184 MPXV genomes of lineage B.1 recovered from GenBank and Nextstrain; and we detected 35 mutations relatively to genome ON563414.3 (a B.1 lineage reference genome). Nonsynonymous mutations occurred in genes encoding central proteins, among which transcription factors and core and envelope proteins, and included two mutations that would truncate a RNA polymerase subunit and a phospholipase d-like protein, suggesting an alternative start codon and gene inactivation, respectively. A large majority (94%) of nucleotide substitutions were G > A or C > U, suggesting the action of human APOBEC3 enzymes. Finally, >1000 reads were identified as from Staphylococcus aureus and Streptococcus pyogenes for 3 and 6 samples, respectively. These findings warrant a close genomic monitoring of MPXV to get a better picture of the genetic micro-evolution and mutational patterns of this virus, and a close clinical monitoring of skin bacterial superinfection in monkeypox patients.

Asfarviruses and Closely Related Giant Viruses.

Acanthamoeba polyphaga mimivirus, so called because of its "mimicking microbe", was discovered in 2003 and was the founding member of the first family of giant viruses isolated from amoeba. These giant viruses, present in various environments, have opened up a previously unexplored field of virology. Since 2003, many other giant viruses have been isolated, founding new families and taxonomical groups. These include a new giant virus which was isolated in 2015, the result of the first co-culture on Vermamoeba vermiformis. This new giant virus was named "Faustovirus". Its closest known relative at that time was African Swine Fever Virus. Pacmanvirus and Kaumoebavirus were subsequently discovered, exhibiting phylogenetic clustering with the two previous viruses and forming a new group with a putative common ancestor. In this study, we aimed to summarise the main features of the members of this group of giant viruses, including Abalone Asfarvirus, African Swine Fever Virus, Faustovirus, Pacmanvirus, and Kaumoebavirus.

Protection from successive Omicron variants with SARS-CoV-2 vaccine and monoclonal antibodies in kidney transplant recipients.

Introduction: Kidney transplant recipients (KTRs) are at high risk of severe COVID-19, even when they are fully vaccinated. Additional booster vaccinations or passive immunization with prophylactic monoclonal antibodies are recommended to increase their protection against severe COVID-19. Methods: Here, we describe the neutralization of SARS-CoV-2 Delta, Omicron BA.1, BA.2, BA.4, and BA.5 variants, firstly by 39 serum samples from vaccinated KTRs exhibiting anti-spike antibody concentrations ≥264 binding antibody units (BAU)/mL and, secondly, by tixagevimab/cilgavimab. Results: No neutralization was observed for 18% of the KTRs, while serum from only 46% of patients could neutralize the five variants. Cross-neutralization of the Delta and Omicron variants occurred for 65-87% of sera samples. The anti-spike antibody concentration correlated with neutralization activity for all the variants. The neutralization titers against the Delta variant were higher in vaccinated KTRs who had previously presented with COVID-19, compared to those KTRs who had only been vaccinated. Breakthrough infections occurred in 39% of the KTRs after the study. Tixagevimab/cilgavimab poorly neutralizes Omicron variants, particularly BA.5, and does not neutralize BQ.1, which is currently the most prevalent strain. Discussion: As a result, sera from seropositive vaccinated KTRs had poor neutralization of the successive Omicron variants. Several Omicron variants are able to escape tixagevimab/cilgavimab.

Correlation of Redondovirus and Entamoeba gingivalis Detections in the Human Oral Cavity Suggests That This Amoeba Is Possibly the Redondovirus Host.

The virome of the human oral cavity and the relationships between viruses and diseases such as periodontitis are scarcely deciphered. Redondoviruses were reported in the human oral cavity in 2019, including in periodontitis patients. Here, we aimed at detecting redondoviruses and at searching for a potential viral host in human saliva. Non-stimulated saliva was collected between December 2020 and June 2021. These samples were tested using real-time PCR regarding the presence of redondovirus and Entamoeba gingivalis DNA. Similarity searches were performed using BLAST against eukaryotic and prokaryotic sequences from GenBank. The redondovirus DNA was detected in 46% of the 28 human saliva samples. In addition, short fragments of redondovirus genomes were detected in silico within Entamoeba sequences. Finally, Entamoeba gingivalis DNA was detected in 46% of the 28 saliva samples, with a strong correlation between redondovirus DNA and E. gingivalis DNA detections, in 93% of the cases. Regarded together, these findings and previous ones strongly support the presence of redondoviruses in the human oral cavity and their association to E. gingivalis as their likely host.

Epidemiological Description and Detection of Antimicrobial Resistance in Various Aquatic Sites in Marseille, France.

Antibiotic resistance is a worldwide public health concern and has been associated with reports of elevated mortality. According to the One Health concept, antibiotic resistance genes are transferrable to organisms, and organisms are shared among humans, animals, and the environment. Consequently, aquatic environments are a possible reservoir of bacteria harboring antibiotic resistance genes. In our study, we screened water and wastewater samples for antibiotic resistance genes by culturing samples on different types of agar media. Then, we performed real-time PCR to detect the presence of genes conferring resistance to beta lactams and colistin, followed by standard PCR and gene sequencing for verification. We mainly isolated Enterobacteriaceae from all samples. In water samples, 36 Gram-negative bacterial strains were isolated and identified. We found three extended-spectrum ß-lactamase (ESBL)-producing bacteria-Escherichia coli and Enterobacter cloacae strains-harboring the CTX-M and TEM groups. In wastewater samples, we isolated 114 Gram-negative bacterial strains, mainly E. coli, Klebsiella pneumoniae, Citrobacter freundii and Proteus mirabilis strains. Forty-two bacterial strains were ESBL-producing bacteria, and they harbored at least one gene belonging to the CTX-M, SHV, and TEM groups. We also detected carbapenem-resistant genes, including NDM, KPC, and OXA-48, in four isolates of E. coli. This short epidemiological study allowed us to identify new antibiotic resistance genes present in bacterial strains isolated from water in Marseille. This type of surveillance shows the importance of tracking bacterial resistance in aquatic environments. IMPORTANCE Antibiotic-resistant bacteria are involved in serious infections in humans. The dissemination of these bacteria in water, which is in close contact with human activities, is a serious problem, especially under the concept of One Health. This study was done to survey and localize the circulation of bacterial strains, along with their antibiotic resistance genes, in the aquatic environment in Marseille, France. The importance of this study is to monitor the frequency of these circulating bacteria by creating and surveying water treatments.

Simple Wastewater Preparation Protocol Applied to Monitor the Emergence of the Omicron 21L/BA.2 Variant by Genome Sequencing.

Detecting and monitoring viruses in wastewater samples have been reported as useful ways of tracking SARS-CoV-2 epidemic trends. However, there is currently no unanimously recognised method of processing samples to identify and quantify SARS-CoV-2 variants in wastewater. We aimed to implement a method that was as simple as possible in order to be used universally. In a study performed between January 2022 and June 2022 in the city of Marseille, France, we first evaluated the impact of the sample preservation strategy. We then compared ultracentrifugation to ultrafiltration and several steps of filtration to determine the optimal approach for virus concentration. As a proof-of-concept, the definitive protocol was applied to next-generation sequencing of SARS-CoV-2 in wastewater to monitor the emergence of the Omicron variant in the city. For sewage water to be processed in the week following the sampling, storage at +4 °C is sufficient, with less than 1 Ct loss. Filtration with a 5 µm syringe filter, then with a 0.8 µm filtration unit, followed by ultrafiltration was the optimal protocol, leading to an average increase of 3.24 Ct when the starting Ct was on average 38 in the wastewater. This made it possible to observe the emergence of the Omicron 21L/BA.2 variant after Omicron 21K/BA.1 by genome sequencing over a period ranging from 20 February to 10 April 2022 in agreement with observations based on patient data. To conclude, by using a simple method requiring only basic filters and a centrifuge as equipment, it is possible to accurately track the relative incidence rates and the emergence of SARS-CoV-2 variants based on sewage samples.

The emergence, spread and vanishing of a French SARS-CoV-2 variant exemplifies the fate of RNA virus epidemics and obeys the Mistigri rule.

The nature and dynamics of mutations associated with the emergence, spread, and vanishing of SARS-CoV-2 variants causing successive waves are complex. We determined the kinetics of the most common French variant ("Marseille-4") for 10 months since its onset in July 2020. Here, we analyzed and classified into subvariants and lineages 7453 genomes obtained by next-generation sequencing. We identified two subvariants, Marseille-4A, which contains 22 different lineages of at least 50 genomes, and Marseille-4B. Their average lifetime was 4.1 ± 1.4 months, during which 4.1 ± 2.6 mutations accumulated. Growth rate was 0.079 ± 0.045, varying from 0.010 to 0.173. Most of the lineages exhibited a bell-shaped distribution. Several beneficial mutations at unpredicted sites initiated a new outbreak, while the accumulation of other mutations resulted in more viral heterogenicity, increased diversity and vanishing of the lineages. Marseille-4B emerged when the other Marseille-4 lineages vanished. Its ORF8 gene was knocked out by a stop codon, as reported in SARS-CoV-2 of mink and in the Alpha variant. This subvariant was associated with increased hospitalization and death rates, suggesting that ORF8 is a nonvirulence gene. We speculate that the observed heterogenicity of a lineage may predict the end of the outbreak.

Matrix-Encoding Gene Diversity of 624 Influenza A/H3N2 Genomes Does Not Show Association with Impaired Viral Detection by Commercialized qPCR Assays.

As for the case of SARS-CoV-2, genome sequencing of influenza viruses is of potential interest to raise and address virological issues. Recently, false-negativity of real-time reverse transcription-PCR (qPCR) assays that detect influenza A/H3N2 virus RNA were reported and associated with two mutations (A37T and C161T) in the Matrix-encoding (M1) gene located on viral segment 7. This triggered a national alert in France. The present study sought to assess the association between the presence of these mutations and potential false negative results of influenza A/H3N2 virus RNA detection by commercialized qPCR assays at the clinical virology laboratory of our university hospitals in southern France. This study focused on the genetic diversity in the M1 gene and segment 7 of 624 influenza A/H3N2 virus genomes obtained from respiratory samples having tested qPCR-positive with M1 gene-targeting assays in our clinical virology laboratory. A total of 585 among the 624 influenza A/H3N2 virus genomes (93.7%) were of clade 3C.2a1b.2a.2, and 39 (6.3%) were of clade 3C.2a1b.1a. M1 gene substitutions A37T and C161T were both present in 582 (93.3%) genomes, only of clade 3C.2a1b.2a.2. Substitution A37T was present in 621 (99.5%) genomes. Substitution C161T was present in 585 genomes (93.8%), all of clade 3C.2a1b.2a.2. Moreover, 21 other nucleotide positions were mutated in ≥90% of the genomes. The present study shows that A37T/C and C161T mutations, and other mutations in the M1 gene and segment 7, were widely present in influenza A/H3N2 virus genomes recovered from respiratory samples diagnosed qPCR-positive with commercialized assays.

Quasispecies Analysis of SARS-CoV-2 of 15 Different Lineages during the First Year of the Pandemic Prompts Scratching under the Surface of Consensus Genome Sequences.

The tremendous majority of SARS-CoV-2 genomic data so far neglected intra-host genetic diversity. Here, we studied SARS-CoV-2 quasispecies based on data generated by next-generation sequencing (NGS) of complete genomes. SARS-CoV-2 raw NGS data had been generated for nasopharyngeal samples collected between March 2020 and February 2021 by the Illumina technology on a MiSeq instrument, without prior PCR amplification. To analyze viral quasispecies, we designed and implemented an in-house Excel file ("QuasiS") that can characterize intra-sample nucleotide diversity along the genomes using data of the mapping of NGS reads. We compared intra-sample genetic diversity and global genetic diversity available from Nextstrain. Hierarchical clustering of all samples based on the intra-sample genetic diversity was performed and visualized with the Morpheus web application. NGS mapping data from 110 SARS-CoV-2-positive respiratory samples characterized by a mean depth of 169 NGS reads/nucleotide position and for which consensus genomes that had been obtained were classified into 15 viral lineages were analyzed. Mean intra-sample nucleotide diversity was 0.21 ± 0.65%, and 5357 positions (17.9%) exhibited significant (>4%) diversity, in ≥2 genomes for 1730 (5.8%) of them. ORF10, spike, and N genes had the highest number of positions exhibiting diversity (0.56%, 0.34%, and 0.24%, respectively). Nine hot spots of intra-sample diversity were identified in the SARS-CoV-2 NSP6, NSP12, ORF8, and N genes. Hierarchical clustering delineated a set of six genomes of different lineages characterized by 920 positions exhibiting intra-sample diversity. In addition, 118 nucleotide positions (0.4%) exhibited diversity at both intra- and inter-patient levels. Overall, the present study illustrates that the SARS-CoV-2 consensus genome sequences are only an incomplete and imperfect representation of the entire viral population infecting a patient, and that quasispecies analysis may allow deciphering more accurately the viral evolutionary pathways.

SARS-CoV-2 quasi-species analysis from patients with persistent nasopharyngeal shedding.

At the time of a new and unprecedented viral pandemic, many questions are being asked about the genomic evolution of SARS-CoV-2 and the emergence of different variants, leading to therapeutic and immune evasion and survival of this genetically highly labile RNA virus. The nasopharyngeal persistence of infectious virus beyond 17 days proves its constant interaction with the human immune system and increases the intra-individual mutational possibilities. We performed a prospective high-throughput sequencing study (ARTIC Nanopore) of SARS-CoV-2 from so-called "persistent" patients, comparing them with a non-persistent population, and analyzing the quasi-species present in a single sample at time t. Global intra-individual variability in persistent patients was found to be higher than in controls (mean 5.3%, Standard deviation 0.9 versus 4.6% SD 0.3, respectively, p < 0.001). In the detailed analysis, we found a greater difference between persistent and non-persistent patients with non-severe COVID 19, and between the two groups infected with clade 20A. Furthermore, we found minority N501Y and P681H mutation clouds in all patients, with no significant differences found both groups. The question of the SARS-CoV-2 viral variants' genesis remains to be further investigated, with the need to prevent new viral propagations and their consequences, and quasi-species analysis could be an important key to watch out.